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Abstract 

Electromagnetically-induced transparency and light storing are studied in 
the case of a medium of atoms in a double A configuration, both in terms of 
dark- and bright-state polatitons and atomic susceptibility. It is proven that 
the medium can be made transparent simultaneously for two pulses following 
their self-adjusting so that a condition for an adiabatic evolution has become 
fulfilled. Analytic formulas are given for the shapes and phases of the trans- 
mitted/stored pulses. The level of transparency can be regulated by adjusting 
the heights and phases of the control fields. 
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I. INTRODUCTION 



It is well known that an atomic medium irradiated by a control laser field may become 
transparent for a signal field which in the absence of the first field would be almost immedi- 
ately absorbed. This phenomenon, known as an electromagnetically-induced transparency 
(EIT) [1,2], has recently been used to drastically change the velocity of a light pulse or even 
to stop or store it: by changing the control field in time it is possible to make the medium 
opaque at the moment at which the signal pulse is inside [3-5]. The pulse is then trans- 
formed into an atomic coherent excitation which is rather robust against relaxation and after 
quite a long time in the atomic scale it is possible to switch the control field on again and to 
release the trapped signal. Such processes have been observed experimentally and explained 
theoretically both in the language of so-called polaritons, being collective atom-|-field exci- 
tations, and in terms of atomic susceptibily [6,7]. An elementary atomic systems for which 
such processes are possible is an atom with three active (resonantly coupled) states in the 
A configuration. 

Adding a fourth active state and a second control field i.e. extending the atomic system 
to a double A configuration allows one to consider new nonlinear, resonantly enhanced opti- 
cal processes [8]. In particular it is possible to simultaneously propagate two optical pulses 
of different frequencies through a medium in the conditions of EIT. In our earlier paper [9] 
we have pointed out that one can stop one pulse and release a pulse of a different frequency 
or two different pulses. We have also mentioned that making the medium transparent si- 
multaneously for two pulses is possible provided that the fields are in some special relation, 
which allows for an adiabatic evolution of the system. If on the other hand the initial condi- 
tions do not satisfy this relation the pulses are partially absorbed, which is connected with 
nonadiabatic phemomena. In the present work we give a detailed quantitative analysis of 
the conditions of a joint medium transparency for two pulses. A discussion in terms of dark- 
and bright-state polaritons, presented in the next section, includes simple analytic formulas 
which allow one to predict the amplitudes, phases and time evolution of the transmitted 



2 



and restored signals and/or the value and space distribution of the atomic coherence due to 
the trapped pulses. Section III shows how the non-adiabaticity of the evolution is reflected 
on the atomic susceptibility. Section IV contains a comparison of the predictions of the 
polariton approach with the results of complete numerical solutions of the Bloch-Maxwell 
equations. 



II. POLARITONS 

We consider a medium composed of atoms in a double A configuration presented in Fig. 1. 
Two driving fields €2,4 (t), induce transparency for two weak signal fields £13. The evolution 
of the atomic density matrix a and the propagation of the signal pulses are described by 
the a set of Bloch-Maxwell equations, while propagation effects for the driving fields are 
neglected. 

In the rotating-wave and slowly-varying-envelope approximations, in the resonance con- 
ditions without relaxations for the transitions shown in Fig. 1, we are left in the first order 
with respect to the signal fields with the equations 

d d 

^^S,^-KliR, + U2a,,), (1) 
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—S3 = -kI{R3 + U^abc), 



—CFbc = U2S1 + UIS3. 



d_ 
dt 



In the above formulas 3 = i^-ih^ ' '^j frequency of the field j, N \s the density of 

atoms, eg is the vacuum electric permittivity, dj are the transition matrix elements: di = dab, 

d2 — dac, o?3 — ddb, di — ddc, Rj — 2!;^) ^2,4 = 2hKi'3 ' ~ —ii^ic^ba, S3 — —in3(^bd and 
a — a{z, t) is the atomic density matrix after transforming-off the rapidly oscillating terms. 

The adiabatic approximation would consist in setting the time derivatives of equal to 
zero, which would mean an evolution during which the atomic upper states are not populated 



at all. However, one can see from the third and fourth equation in the set (1) that this is 
possible if the four laser fields remain in a certain proportion: = This reflects the 
fact that during an adiabatic evolution each atom makes a continuous transition from the 
initial state b to the dark state, which is a superposition of the two lower states b and c with 
the coefficients determined by the instantaneous values of the intensities of the signal and 
control fields. However, in a double A system this combination is a dark state simultaneously 
for the two A's only in the special situation. 

After eliminating the variables 81^3 the above equations can be rewritten in a new set of 
variables - the so-called polaritons 

Ri — exp(^ arg(t/2)) (cos 6 cos + sin 6 cos 0$ + sin (j)X) , 

i?3 = exp(iarg(C/4))(cos^sin0^ + sin ^ sin 0$ — cos(/>X), (2) 

(Jbc = - sin 6'^' + cos 6*$, 

where sin^ = (It/sP + W + 1)"^^^ and tan0 = 

The evolution equations for the polaritons in the case of time-independent driving fields 
read 

d d d 

—'^ + ccos^9—'i> + csm9cos9—^ = 0, 

ot oz oz 

- sin^ 6 d"^ d"^ 1 

Z-^* + = — ^{i^l cos^ <^ + ^^3 sin^ <t>)^ - sin cos - kI)X, (3) 

cosc^ or oV smfc' 

9 d sine/) cos 2 / 2 • 2 , , 2 2 , x ^ 

In the case of a single A system (i.e. for X = 0) an adiabatic evolution meant that the dark 
state polariton traveled with the velocity ccos^^ keeping its shape, with the bright state 
polariton $ = 0. The only necessary condition of adiabaticity was that the pulses 1 and 
2 should be smooth enough. For a double A system, again assuming smoothness of all the 
pulses, the new element is that the bright state polariton X is in general different from zero 
at the beginning of propagation: the decomposition of the incoming signal fields is given by 
Eqs (2) (with $ = 0) 
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Ri — exp[i arg(t/2)o] (cos 9o cos (/>o^ + sin (l)oX), 

i?3 = exp[zarg([/4)o](cos6'osin(/)o^' — cos0o^)- (4) 

Then during the evolution ^ keeps its shape while X is damped. Thus the conditions of 
adiabaticity are gradually created and the three dynamical variables tend to the following 
values 

Ri exp[i arg(?72)] ^ ^"'^ ^ (cos 0o exp[-i arg([/2)o]it!i + sin 0o exp[-i arg(?74)o]i?3) 

cos UQ 

cos sm (p 

R3 exp[zarg(t/4)] — (cos0oexp[-iarg(f/2)o]i?? + sin ^oexp[-z arg([/4)o]i?3), (5) 

cos Uq 

sm 

CTbc — (cos0oexp[-iarg(C/2)o]-R? + sin0oexp[-iarg(C/4)o]i?3). 

cos Oq 

The argument of the functions i?? 3 in the r.h.s of Eqs (5) has to take into account the shift 
by /o ccos^ 6{T)dT during the time t. In particular if we assume that at t = the edge of the 
pulse reaches the sample, when calculating the position of the maximum we have to take 
into account that the maximum moves first with the velocity c until it reaches the sample 
and later travels with the velocity ccos'^ 6. 

The above formulas allow one to predict the shape and the numerical parameters of the 
transmitted pulses in typical EIT, when the control fields are kept constant, as well as in 
the process of light storing, when the control pulses are slowly varying. Note in particular 
that a single initial signal in the presence of two control fields leads to an appearance of the 
other signal pulse. In Section IV we will compare the predictions based on these formulas 
with the results of the numerical solutions of the Bloch-Maxwell equations. 

III. ATOMIC SUSCEPTIBILITY 

A complementary picture of the above evolution is obtained in terms of the susceptibility. 
If we rewrite the first three equations of Eq. (1) in the original variables admitting the 
detuning and relaxations terms 



ma = ^^3^ + ^4(^6c + AgfTbd (6) 

where ^2,4 = ^^if^, ^Ai = - - hui - ir«, hA^ = E^ - E^ - tiw^ - fF'^ and 
hS — Eh + hcol — Ec — hco2 — are the detunings including the relaxation rates for the 
coherences ab, db and be. 

If we now pass to the frequency domain, assuming that u>i — UJ2 — UJ3 — uj^, we can 
calculate the elements of the density density matrix a and express the components of the 
polarization in terms of the signal fields 

Ndiabaii^) = 27reo[xii(w)ei(w) + Xi3(c^)e3(t^)] 

Ndsabdiuj) = 2-Keo[x3i{uj)ei{u;) + X33{i^)^3{(^)], (7) 

where 
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Xsil^J- 47r^e^^ (Ai-cu)(A3-cc;)(<5-tj)-(Ai-a;)|fi4P-(A3-a;)|02P' ^ 

X33(^) = 

It is important to notice that under the assumptions of the polariton analysis, i.e. in the 
resonance conditions and neglecting the relaxations, the expressions for the susceptibility x 
become singular at a; = 0, e.g., 

Xii(^) = 

_ Nd.di n2ni 

^''^'^^ 47:heoUJ^-uj{\n2\^+\n4^y 
with analogous expressions for X3i and xas- One cannot thus speak of a transparency window 
and the pulses' propagation occurs with a significant distortion. Fulfilling the adiabaticity 
condition ^ — ^ (cf. the third and fourth of Eqs (1)) means that the singularities in Eqs 
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(7) cancel out and the induced transparency inside a transparency window of a finite size is 
possible. In that case after substituting 63 = ^i^i^jf we obtain 

Ndiaba^'2n€oXii{io)€i{u), (10) 

with 

^''^""^ - Ai^he, a;2-(|02|2 + |04|2)' ^^^^ 

which is the expression as for a single A system, with the only difference that the denomi- 
nator is corrected by the |r24|^term, which means that the transparency window is widened 
compared with the case of a single A system. 



IV. NUMERICAL ILLUSTRATION 

A numerical illustration of the above results is presented below. We have performed 
calculations for the atomic model with £'a=-0.10 a.u., £'6=-0.20 a.u, £'c=-0.18 a.u., E^—- 
0.05 a.u., the width of the upper levels a and d connected with the spontaneous transitions 
to the lower levels h and c were taken V^l = 2 Ax 10"^ a.u., from which the dipole moments 
have been calculated; the relaxation of the lower-states' coherence has been neglected. The 
medium density has been taken N = 3x 10"^^ a.u. (2 x 10^^ cm~^) . The length of the sample 
was 10^ a.u (0.5 mm). The initial length of the signal pulses was 10^^ a.u. (2.4 ^s) and their 
amplitudes, modeled by a sine-square, were of order of 10~^° a.u., which correspond to the 
power density of 3.5 x 10~^ Wcm""^. The amplitudes of the control pulses were larger by 
about an order of magnitude. 

Because there are essentially two reasons of nonadiabaticity of the evolution, namely 
the discontinuity of pulses or the fields' faihng to satisfy the proportion following Eq. (1), 
we first check the role of the former effect. In Fig. 2 we show the shape of the initially 
rectangular pulse in a single A system in a control field with a time-independent amplitude. 
The evolution is initially nonadiabatic. However, the signal propagating in the medium is 
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gradually smoothed, which corresponds to an absorption of the polariton An analogous 
calculation with a smooth (sine square) pulse yields to a very good approximation a conser- 
vation of the pulse shape during the propagation, which means an absence of the polariton 
$ from the very beginning. Because in our further investigations on double A systems we 
use sine square pulses, it follows that we may indeed neglect the polariton as mentioned 
in Section II. 

In Fig. 3 we show the transmitted pulses 1 and 3 compared with their initial shapes. 
The horizontal lines show the values calculated from Eq. (5). The predictions of the heights 
of the pulses are excellent. We have checked that Eq. (5) predicts properly also the phases 
of the transmitted signal pulses. 

The effective transparency for the two pulses can be regulated by choosing the heights 
and phases of the control fields. In Fig. 4 we present the transmitted pulses 1 and 3 for two 
different phase relations. By a sutiable choice of the phase of any of the control fields one 
can regulate the heights of the signal pulses. 

If we switch the control fields off the medium becomes opaque for the two signal pulses, 
which are then "stopped" in the form of the single atomic coherence abc, the form of which 
is given by the third of Eqs (5). Of course the absolute value of this coherence depends 
on the place inside the sample but the phase relations are more general. In Fig. 5 we show 
the argument of the complex coherence as a function of the phase shift of the field 3; the 
simultaneous adiabatic switch-off of the two control fields has been modeled by a hyperbolic 
tangent. One can see that the predictions of the last of Eqs (5) are also excellent. 

In Fig. 6 we show the space distribution of the coherence due to the trapped pulse. 
The length of the sample was 3.5 x 10^ a.u. and the control pulses were switched off again 
simultaneously as a hyperbolic tangent, but somewhat later than before so that the whole 
pulse could be trapped inside the medium without any additional effects of the boundary 
of the sample. The distribution can be modeled with a good accuracy by a sine square, i.e. 
the common shape of the incoming pulses. No parameters have been fitted: the amplitude 
has been obtained from the last of Eqs (5), the width is the original width multiplied by the 
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compression factor cos^ (the pulse has been compressed when entering the sample but did 
not change its width any more during the storing stage) and the position of the maximum 
has been calculated as explained above following Eqs (5). 

V. CONCLUSIONS 

We have discussed a simultaneous propagation of two signal pulses in the medium of 
four-level atoms in the double A configuration. We have shown quantitatively that it is 
practically possible to divide the incoming pulses into a dark- and bright-state polaritons, 
of which the dark one survives during the evolution. This is connected with the pulses' self- 
adjusting during the initial stage of the propagation, so that the condition of adiabaticity 
becomes fulfilled. A knowledge of the shape and position of the dark-state polariton allows 
one to predict the characteristics of the transmitted pulses or, in case the pulses were stored, 
the value and space distribution of the induced atomic coherence. Those characteristics 
can be dynamically controlled by changing the parameters of the control fields, e.g. their 
phase difference. We have also demonstrated that an adiabatic character of the the pulses' 
evolution is reflected in a cancellation of singularities of the atomic susceptibility. 
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FIGURES 

FIG. 1. Level and coupling schemes; the indices 1 and 3 refer to signal fields and 2 and 4 - to 

control fields. 

FIG. 2. The shape of the propagating initially rectangular pulse as a function of the local time 
t' = t — ^, in the case of a single A system for a control field of a time- independent amplitude 
€2 = 1.2 X 10~^ a.u.: at the beginning of the sample (1), at the depths: 3 x 10^ a.u. (2), 1.5 x lO''' 
a.u. (3), 3 X 10^ a.u. (4). Note the process of smoothing the pulse. 

FIG. 3. The shapes of the transmitted pulses 1 (curve 1) and 3 (curve 3), compared with 
their initial shapes (curves 10 and 30, respectively, for the control fields €2 = 1.2 x 10~^ a.u. and 
e4 = 1.8 X 10~^ a.u. The horizontal lines show the predicted heights of the transmitted pulses. 
Note also the slowdown of the transmitted pulses. 

FIG. 4. The shapes of the transmitted pulses 1 (curves la and lb and 3 (curves 3a and 36) for 

the data as in Fig. 3 but with additional phase shifts (pj of the pulse j: ^!'2 = f , 03 = ^4 = 
(case o) and (?!>4 = | (case b). 

FIG. 5. The phase shift of the coherence due to the trapped pulses as a function of the 
phase shift of the field 3. The crosses are the results of the numerical computations and the solid 
line is obtained from the last of Eqs (5). The control fields of the maximum values as in Fig. 3 
were switched off simultaneously as a hyperbolic tangent. 

FIG. 6. The spatial distribution of the coherence due to the trapped pulse. The crosses are 
the results of the numerical computations and the solid line is a sine square, of which the amplitude 
has been calculated from the last of Eqs (5), the position of the maximum is given by the formula 
following Eqs (5) and the width is equal to the width of the incoming pulses divided by cos^ 60 . 
The control fields were of the same shape as in Fig. 5 but were switched off somewhat later. 
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